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We investigate the relationship between colour confinement and the monopoles derived from the 
Cho-Duan-Ge decomposition. These monopoles, unlike Dirac and 't Hooft monopoles, do not 
require a singular gauge field and are defined for any choice of gauge (and are not just restricted 
to, for example, the maximum Abelian gauge). The Abelian decomposition is defined in terms of 
a colour field n; the principle novelty of our study is that we have used a unique definition of this 
field in terms of the eigenvectors of the Wilson Loop. This allows us to investigate the relationship 
between the gauge invariant monopoles and confinement both analytically and numerically, as 
well as retaining the maximal possible symmetry within the colour field so that it is able to see all 
the monopoles in an SU(Nc) calculation. 

We describe how the Abelian decomposition is related to the Wilson Loop, so that the string ten- 
sion may be calculated from the field strength related to the decomposed (or restricted) Abelian 
field. We demonstrate that for an area law scaling of the Wilson Loop there must be discontinu- 
ities in the restricted field, and discuss the structures in the colour field which may cause these 
discontinuities, which turn out to be magnetic monopoles. If these monopoles are present, they 
will lead to an area law scaling of the Wilson Loop and thus be at least partially responsible for 
confinement. 

We search for these monopoles in quenched lattice QCD. We show that the string tension is 
dominated by peaks in the restricted field strength, at least some of which are located close to 
structures in the colour field consistent with with theoretical expectations for the monopoles. We 
show that the string tension extracted from the monopole contribution to the restricted field is 
close to that of the entire original field; again suggesting that confinement can at least partially be 
explained in terms of these monopoles. 
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1. Introduction 

The mechanism behind quark confinement is an enduring problem in QCD. Although sev- 
eral models have been proposed - for example, center vortices and a dual Meissner effect 
due to magnetic monopoles ^ - none have yet been demonstrated convincingly. Here we study 
the possibility of confinement due to monopoles constructed from the Cho-Duan-Ge (CDG) de- 
composition (sometimes referred to as the Cho-Faddeev-Niemi decomposition) [Q]. Unlike Dirac 
monopoles and 't Hooft (Maximum Abelian Gauge) monopoles, the CDG decomposition respects 
the gauge symmetry and does not require a singular gauge field. The decomposition is constructed 
from a colour field, n, which is built from a matrix 6. Recent work [Q] has demonstrated that 
the monopole contribution dominates the confining string, using monopoles constructed from one 
particular choice of 6 G SU (Nc) /U (Nc — 1); however in this case only one of the possible Nc — 1 
types of monopole is visible. Here we consider a different choice of d G SU(N c )/(U(l)) Nc -\ and 
investigate whether this the monopoles apparent in this construction may also lead to confinement. 
Our full results and methods will be published later In section ^ we discuss the Abelian decom- 
position and its relation to the Wilson Loop and thus quark potential; in section || we discuss how 
monopoles may arise in this construction and lead to confinement; we present numerical evidence 
in section ^| and conclude in section [|. 

2. Abelian decomposition and Stokes' theorem 

The confining potential in an SU(iVc) gauge theory can be measured using the Wilson Loop, 

W L [C s ]~tx(W[C s }) W[C s ]=P[e-^ dx ^] (2.1) 

Nc 

for a closed curve C s of length L which starts and finishes at a position s, where P represents path 
ordering and the gauge field, A^, can be written in terms of the Gell-Mann matrices, X", as A"X a . 
The expectation value of the Wilson Loop scales as (Wl[C s ]) ~ e~~ , where s# is the area of the 
surface enclosed by the curve C s and s is the string tension. We only consider planar Wilson Loops: 
C s is a rectangle of temporal extent T and spatial extent R. The quark-quark potential is given by 
V(R) = \im T ^\og((W L [C s \))/T. 

To circumvent the path ordering, we split C s into infinitesimal segments of length 8a, and 
define the gauge link as U a G SU(N C ) = P[e~ ig ^ +5 ° A ° da ] ~ e - igS ° Aa . < a < L represents the 
position along the curve and we write A a = A /i ( ff ^(jc(o')). We have assumed that the gauge field is 
differentiable. W [C s ] can then be written as 

L-Sa 

W[C S ] = lim U a . (2.2) 

5<T ^°(T=0,5<T,25<7 ) ... 

We proceed by inserting an identity operator between each pair of gauge links along the curve. 
The goal is to replace U with an Abelian field, and evaluate the Wilson Loop using Stokes' theorem. 
Previous work [@] has used the identity / = / d@& \e Nc ) (e Nc \ @\ for G SU(N C ) /U(N C - 1) and 
€pj c a unit colour vector; but we choose differently, with no need for the integral and with in a 
larger group so that it can 'see' more types of monopole. 
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We introduce a field d a = d(x(a)), which, for the moment, we shall take to be an element 
of U(Nc), at each point along C s and insert the identity operator d a dl between each of the gauge 
links. 6 is chosen so that dlu a d a+ s a is diagonal. 6 S therefore contains the eigenvectors of W[Q]: 
W [C V ]0 S - = 6 S D, where D is some diagonal element of SU (Nc). As the phases of the eigenvectors 
are arbitrary, this definition only determines 6 up to a U (Nc) transformation 6 —> 6%. Fixing 
these phases and the ordering of the eigenvalues by some arbitrary fixing condition gives a unique 
choice of 6 € SU(N c )/(U(l)) Nc ~ 1 . Under a gauge transformation U a — > A a U a A^ Sa for A G 
SU (N c ), d ^Ad%, where the U (N c ) factor % depends on the fixing condition. With Blu a B a +8a = 

g'^A-/ diagonal 

d}W [C s ] 6 S = e iL v *»»*> v fcs d W ; (2.3) 

removing the non-Abelian structure and the path ordering. 

We may extend this definition of 6 across all space by constructing nested curves in the same 
plane as C s and then stacking these curves on top of each other in the other dimensions. We then 
define 6 so it diagonalises W along each of these curves. We now introduce a restricted S\J(Nc) 
gauge field Up which is diagonalised by d, so, for A 7 a diagonal Gell-Mann matrix 

[X J \ dxUnjdx+jiSo] =0 

%A+Soffil*-»i =0 < =W0l (2.4) 

is satisfied across all of space-time. We can introduce a second field X such that U^(x) = XnUn, 
and restrict X^ by imposing the condition 

tr[ni(Xl x -X^ x )}=0 V diagonal. (2.5) 



If there are multiple solutions to equations ( [2.4] ) and ( |2.5[ ) we select the solution which maximises 
tr(X). Under a gauge transformation, n x — > A x n x A\, (x) — > A x U^ ltX A^ x ^g a and X^^ — > A x X^ x Al, 
so equations (2.4) and (^5f) are gauge invariant. Equations (2.4) and (2.5) are the lattice versions 
of the defining equations of the CDG decomposition^]], which in the continuum is described by 

A fl =A fl +X fl D M [AK'=0 
=tr(n j X) [A]a =d^a - ig[A, a] 



Xj diagonal 



n J tr(n j A IJL ) - -ig l [n J ,d^n J ] 



(2.6) 



The corresponding field strength is f^ v [A] = n^F^y with the gauge invariant 

F^y = j \dfitt(n^A v ) — d v tr(n->An)] — ^tr(nj[d^n k ,d v n k ]). We express the restricted field as = 

^ ^ since £ = jj ^ cmye ^ wg ^ ^ fl^J/] =W[C S ,U] = 

8 S W[C S , dWd]dJ. Applying Stokes' theorem to the Abelian field d^U^ rX d x+jlSa gives 

e}W[C s ]d s = e ! '^ J ^nal^7^MV^V ) (2 .7) 

where dS^ v is an element of the planar surface bounded by C. Thus we may expect structures in 
the field strength Fj{ v to lead to a confining potential. 

Whenever U is differentiable, F^ v is an exact derivative: Fj{ v = d^(u J v ) — d v (u 3 ^) An area 
law scaling of the Wilson loop requires that U (and thus 6) is non-differentiable at certain points. 
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Lattice size (lattice units) 


Spatial Lattice size (fm) 


P 


lattice spacing (fm) 


Number of configurations 


I6 J 32 


2.30 


8.0 


0.144(1) 


91 


I6 3 32 


1.84 


8.3 


0.115(1) 


91 


16 3 32 


1.58 


8.52 


0.099(1) 


82 


20 3 40 


2.30 


8.3 


0.115(1) 


20 



Table 1: Parameters for our simulations. 



3. Parametrisation and divergences of the field 

d s is denned as the matrix of eigenvectors of W [C s ] . It is discontinuous when W [C s ] is dis- 
continuous or has degenerate eigenvalues, and in another situation: in SU(2), we parametrise 6 
as 

e=( cosa . < Wc W, (3-D 

\ isinae lc cos a I 

with c £ R, < a < n/2 and dj determined by the fixing condition. At both a = and a = n/2, 
the parameter c is ill-defined, and points where a = n/2 can lead to monopoles. The parameter c 
may wind itself around these points, creating a discontinuity in B. In the plane of the Wilson Loop, 
we parametrise space-time in polar coordinates (r, 0), with the origin at the point where a = n/2. 
At some infinitesimal radius r, continuity of 6 demands that c(r, 0=0)= c(r, (j) = 2n) + 2nv for 
integer winding number v, and if c is ill-defined at r = we may find v which implies a 
discontinuity in 6. We may evaluate F^ v using an integral of u around a loop of infinitesimal r 
centred at r = 0: F^ v = 8(x) § d&u J s , and it can be shown that % is proportional to dgc This 
gives a 5-function in Fj{ v : a CFN magnetic monopole. It is reasonable to expect that the number of 
these monopoles will be proportional to the area of the Wilson Loop, leading to an area law scaling. 

In SU(3), we parametrise 6 in terms of six variables < a\ < n/2, < ai < n/2, < < 
n/2, c\ G R, C2 £ R, and C3 £ R, and two fixed constants, d$ and d&, as 



\ / cos t/2 ; sin«2^ ,c - \ / cosaj 1 sin^j e ,L I \ • / 1 3 ■ j8 1 8 

d = [ cosn, ,sin a3f "3 1 ,sm fllf --l cos», ) e 3 ^ +U ' A . (3.2) 

isintf^e - "^ cosfl3 / V isina 2 c JC 2 cos« 2 / V 1 



There will be monopoles when one of the a,- is n/2 and one of the c\ winds itself around that point. 

The area law scaling of (tr(Wj,[C])) may thus be related to 5-functions in the reduced field 
strength caused by discontinuities in the 0-field at a = n/2 (seen in the magnetic, tr(n[^n,5 v «)]> 
term within F^ v ) and characterised by a non-zero winding of the parameter c around these monopoles. 



4. Numerical results 

We generated 16 3 32 and 20 3 40 quenched configurations in SU(3) with a Tadpole Improved 
Luscher-Weisz gauge action [^] using a Hybrid Monte Carlo routine [[7|] (see table |]). The lattice 
spacing was measured using the string tension a ~ (420MeV) 2 . We fixed to the Landau gauge and 
applied ten steps of improved stout smearing [Q] with parameters p = 0.015 and £ = 0. 6 and U 
were calculated numerically. 

In figure |T] (top plot), we plot a slice of the restricted field strength at fixed y and z coordinates, 
and see that it is indeed dominated by peaks a few lattice spacings across. Plots on neighbouring 
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Figure 1: A cross section of the field strength F^ at z = 0,y = 10 on one configuration (top), and the U string 
tension excluding peaks of height |F xt | > c gp^ max ^ from the average over Wilson loops (bottom). Due to the 
limited resolution of the lattice, the extrapolated contour lines have an error of up to one lattice spacing. 



slices of the lattice show a very different pattern, suggesting that these are point-like objects rather 
than strings or membranes. In figure [j] (bottom), we show that these peaks are responsible for the 
string tension by plotting the quark potential only averaging over those Wilson loops in the xt plane 
which exclude these peaks: only including loops where |Fjcf(jc)| < c <oF^ nax \x), Yx within the curve 
C s . F^ max ^ is the configuration maximum value of the field strength. The string tension decreases 
as more of these objects are excluded, suggesting that these maxima indeed cause the area law 
scaling. 

In figure || (top), we investigate whether the peaks in F xt are close to the maxima of a,- (on 
the lattice, these will not be precisely at a = Tt/2). We have only shown data for one of the £?,• 
parameters; those peaks which are not close to a maximum of a in this plot are near the maximum 
of a different a*. The contours lines show the maxima of F#, combining the two components of 
the gauge field, while the level of shading shows the values of au with a dark colour indicating 
at ~ and a light colour a\ ~ 7t/2. The overall picture is a little ambiguous, but consistent with 
expectations. We also show a similar plot for c; (figure || bottom). This should wind itself around 
the maxima, shown by a gradual darkening of the background as we rotate around the peak. Again, 
we only show data for one of the three c, but we see some examples which suggest that c indeed 
winds itself around the peak (for example at the center of the plot). 

In figure || we compare the string tension for the original gauge field U, the decomposed gauge 
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Figure 2: A comparison between the peaks in F^' 8 (red contours) against the peaks in 2a\/n (shaded 
background, top) and the angle C2 (shaded background, bottom) 

field U, and the monopole contribution. The monopole string tension is extracted from the CDG 
decomposition of d^ x = O^Uu^x+au. (rather than the non-gauge invariant 6%6 x+a n which gives 
too noisy data), with L^, Y the gauge field subjected to 600 sweeps of stout-smearing at p =0.1: 
enough smearing will destroy any structure, and any signal will be from the 6 field calculated 
without smearing. The U field is identically equal to the U field on those Wilson Loops used to 
define 6; but to increase statistics, we averaged U over every Wilson Loop; so our results for the U 
and U fields differ. Both the restricted field and the monopole field dominate the string tension. 

5. Conclusions 

We have investigated whether the confining quark potential in quenched QCD is caused by 
CDG gauge-invariant monopoles. Our main novelty is to construct the CDG colour field ?i ; = 
QXiQ* from the eigenvectors of the Wilson loop, which allows us to access the full symmetry 
group of the monopoles and permits a theoretical discussion. Our numerical results suggest that 
the restricted field strength is dominated by peaks one lattice spacing across and that these peaks 
are responsible for the confining potential. We see some, though not yet convincing, evidence 
supporting the theoretical expectation that these peaks are close to the point where the SU(2) sub- 
components of 6 are off-diagonal and that the 6 field winds around these peaks. This suggests that 
these monopoles are at least partially responsible for confinement. We will expand this argument 
and give full details of the calculation in a subsequent publication. 
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Figure 3: The string tension extrapolated to infinite time for the original gauge field U, the restricted gauge 
field U, the over-smeared field U and the monopole field M (left); the fit results for the string tension (right) 
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